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ABSTRACT 

We report the results of a pilot survey for X-ray emission from a newly discovered class of AGB 
stars with far-ultraviolet excesses (fuvAGB stars) using XMM-Newton and Chandra. We detected 
X-ray emission in 3 of 6 fuvAGB stars observed - the X-ray fluxes are found to vary in a stochastic 
or quasi-periodic manner on roughly hour-long times-scales, and simultaneous UV observations using 
the Optical Monitor on XMM for these sources show similar variations in the UV flux. These data, 
together with previous studies, show that X-ray emission is found only in fuvAGB stars. From 
modeling the spectra, we find that the observed X-ray luminosities are ^ (0.002 — 0.2) L e , and the 
X-ray emitting plasma temperatures are ~ (35 — 160) x 10 6 K. The high X-ray temperatures argue 
against the emission arising in stellar coronae, or directly in an accretion shock, unless it occurs on 
a WD companion. However, none of the detected objects is a known WD-symbiotic star, suggesting 
that if WD companions are present, they are relatively cool (< 20, 000 K). In addition, the high 
X-ray luminosities specifically argue against emission originating in the coronae of main-sequence 
companions. We discuss several models for the X-ray emission and its variability and find that the 
most likely scenario for the origin of the X-ray (and FUV) emission involves accretion activity around 
a companion star, with confinement by strong magnetic fields associated with the companion and/or 
an accretion disk around it. 

Subject headings: binaries: general — binaries: symbiotic —stars: AGB and post-AGB — stars: 

mass-loss — stars: individual (EYHya, Y Gem, CIHyi) — circumstellar matter 


1. INTRODUCTION 

Almost all of our current understanding of the late 
evolutionary stages of (^ 1 — 8) M 0 stars is based on 
single-star models. However, binarity can drastically af¬ 
fect late stellar evolution by (a) cutting short normal 
AGB evolution or even preventing stars from reaching 
the AGB, due to a phase of strong binary interaction, 
when the primary was an RGB star, and (b) prolonging 
post-AGB evolution due to mass-transfer back onto the 
primary from a circumstellar disk around the compan¬ 
ion (van Winckel 2003). It has long been argued that 
binarity is responsible, directly or indirectly, for the dra¬ 
matic and poorly-understood changes in the history and 
geometry of mass loss that occurs in stars as they evolve 
off the AGB to become PNs. A variety of binary mod¬ 
els (e.g., review by Balick & Frank 2000) have been pro¬ 
posed, which can lead to the generation of accretion disks 
and magnetic fields. The latter are likely the underlying 
physical cause for the highly collimated jets that have 
been proposed as the primary agents for the formation 
of bipolar and multipolar PNs (Sahai & Trauger 1998, 
Sahai et al. 2011a). 

However, observational evidence of binarity in AGB 
stars is sorely lacking simply because AGB stars are very 
luminous and variable, invalidating standard techniques 
for binary detection such as radial-velocity and photo¬ 
metric variations due to a companion star. Sahai et 
al. (2008: Setal08) therefore used an innovative tech¬ 
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nique of searching for UV emission from AGB stars with 
GALEX (Morrissey et al. 2007) that exploits the fa¬ 
vorable secondary-to-primary photospheric flux contrast 
ratios reached in the UV for companions of spectral type 
hotter than about GO (T e ff=6 000 K) and luminosity, 
L > 1I/Q. Setal08 detected emission from 9/21 objects 
in the GALEX FUV (1344- 1786 A, effective A=1516 A) 
and NUV (1771 — 2831 A, effective A=2267A) bands; 
since these objects (hereafter fuvAGB stars) also showed 
significant UV variability, Setal08 concluded that the UV 
source was unlikely to be solely a companion’s photo¬ 
sphere, and was dominated by emission from variable 
accretion activity. 

From a subsequent search of the MAST0/GALEX 
archive, we found about 100 fuvAGB stars > 5cr de¬ 
tections in the FUV band. We required the UV source 
position to be coincident, within 3", with the optical po¬ 
sition of the AGB target stars (sp. types M4 or later). 
The chance coincidence probability with random sources 
for our objects is extremely low (Setal08). Even for the 
hottest sources in our catalog (sp. type M4), the detected 
FUV fluxes (->20 /iJy) correspond to a significant excess 
above photospheric emission - e.g., for an M4III giant 
with L ~ 6000 Lq, T e ff ^ 3560K (e.g., Ridgway et al. 
1980, Perrin et al. 1993), and distance 500 pc, the black- 
body flux at 1500 A is 6.2 /iJy. Many fuvAGB stars show 
extreme UV variability as well, such as Y Gem (Sahai et 
al. 2011b: Setalll). 

We report here the results of a pilot search for X-ray 
emission from a small subset of fuvAGB stars. Our detec¬ 
tion of X-ray emission from 3 objects more than doubles 
the known number of X-ray emitting AGB stars, pro- 
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viding, for the first time, high-quality X-ray spectra and 
X-ray and UV light curves. 

2. OBSERVATION & RESULTS 

We observed 5 sources with XMM-Newton using 
EPIC in full window mode, (AO-13 Priority C proposal 
072034), and 4 sources with the Chandra X-Ray Observa¬ 
tory (CXO) using ACIS-S (Cycle 15 proposal 15200476.) 
Three sources, CIHyi, Y Gem & EYHya, were common 
to both programs, and were detected in X-ray emission 
(Fig.QJgb). Simultaneous UV observations, providing us 
UV light-curves for these 3 sources were obtained us¬ 
ing the Optical Monitor on XMM through the UVM2 
(A = 231 nm) and UVW2 (A =212 nm) filters in Fast 
Mode. Essential observing details are given in Table[lJ 
Data were reduced following standard procedures, and 
cleaned by removing time intervals affected by high back¬ 
ground, and the cleaned observations were used to ex¬ 
tract the X-ray spectral energy distributions in the 0.3- 
10 keVrange (see, e.g., Sanz-Forcada et al. 2011). 

Short-term flux X-ray variations are seen during the 
course of each observation, spanning a factor ~ (2 — 3.5) 
on hour-long timescales (Fig.[2ji,b). The energy ranges 
used to extract these curves was 0.3 — lOkeV for the 
XMM/EPIC instruments, and 0.3 —7keV in CXO/ACIS 
(for the latter, the removal of data for E > 7 keV is be¬ 
cause that region has more problems due to background 
noise, as explained in Chandra documentation.) For 
Cl Hyi and EY Hya, the variations appear to be stochas¬ 
tic or quasi-periodic. In Y Gem, the light-curves give 
stronger evidence of a discrete period of about 1.2-1.4 hr, 
best seen in the CXO/ACIS-S light curve (Fig.Ebl) that 
spans a longer period than the EPIC one. For the longer 
time-scales between the XMM and CXO observational 
epochs for each source, we find changes in the observed 
fluxes up to a factor 1.8 (in YGem and EYHya); how¬ 
ever a direct comparison cannot be made between the 
spectra as these were taken with different instruments. 

The UV fluxes also show short-term variations like 
the X-ray fluxes (Fig. [2^). Since the UVM2 and UVW2 
data were taken sequentially (the OM can take data only 
through one filter at a time), we have scaled the UVW2 
count rate to match that of UVM2 at the time of tran¬ 
sition between the two. We find that the UV and X-ray 
variations appear reasonably well correlated in the case 
of Y Gem and Cl Hyi. For EY Hya, the correlation ap¬ 
pears to be weaker. 

The ISIS package (Houck & Denicola 2000) and the As¬ 
trophysics Plasma Emission Code (APEC, Smith et al. 
2001) were used to fit the spectra with thermal models. 
The background spectrum was estimated from an empty 
field on the same detector chip as the source, and pro¬ 
vided to the ISIS software that then used it to compute 
and fit a model ” Source+Background” spectrum to the 
total spectrum extracted from a circular aperture cen¬ 
tered on the source. For XMM, we simultaneously fit 
the data from the EPIC pn, MOS 1 and MOS 2 detec¬ 
tors. For each object we first fit the EPIC spectrum as 
it has higher S/N than the ACIS one. The elemental 
abundances used are those of Anders & Grevesse (1989), 
except for the Fe abundance ([Fe/H]=logl0 of the Fe 
abundance divided by the solar value) which was chosen 
to be a free parameter only if the S/N was sufficiently 
high, or fixed appropriately (see below). 


All 3 sources also show a weak line feature at ^ 
6.3 — 6.8 keV in their XMM spectra, that is likely due to 
emission from (a) the Fe XXV, XXVI lines complexes at 
~ 6.7keV from the very hot plasma that we have found 
in these objects, and/or (b) the 6.4 keV Fel Ka line. 
We present an expanded version of the Y Gem XMM pn 
spectrum that clearly shows these lines in Figure[3l The 
origin of Fe I Ka line may be the same as in Young Stellar 
Objects (YSOs), where it has been inferred to be floures- 
cent emission from cold, neutral (< 1 MK) material in a 
disk irradiated with energetic (E > 7.1 keV) X-rays (Fa- 
vat a 2005). 

The most robust parameters from our models are 
the X-ray temperature (Tx) and the emission measure 
(EM), followed by the hydrogen absorption column den¬ 
sity (NH) - these provide us with the intrinsic X-ray 
luminosity (Lx) of each of our objects. We find that 
Lx - (0.002 - 0.2) Lq and Tx - (35 - 160) MK for our 
detected objects (Table[2]). For the minor Fel line com¬ 
ponent, our derived values of the line flux are sensitive 
to the emission measure (EM) of the hot plasma that 
contributes to the nearby Fe XXV, XXVI complex, and 
are less robust than indicated by the formal uncertainties 
derived from the least-squares fitting. 

The simplest fit to the Y Gem XMM spectrum (Epoch 
2, Fig. Util), which has the highest S/N of all the spectra 
reported here, requires two APEC components with the 
same temperature but with different levels of absorption 
(NH(1) and NH(2)) and emission measure (EM(1) and 
EM(2)), and a gaussian Fel line. For the latter, we had 
to impose the condition that the line is narrow, and we 
therefore fixed its width (FWHM) at a nominal value 
of 0.002 keV0; the line energy was fixed at 6.4 keV. We 
allow [Fe/H] to vary. The Fel line is assumed to be 
extincted by NH(2), which is significantly smaller than 
NH(1) - assuming this line comes from a cold disk, the 
lower attentuation could be an effect of disk inclination 
and viewing geometry, with the flourescing region on the 
disk having a smaller amount of extincting gas between it 
and the observer. The derived [Fe/H] is not significantly 
different from the solar value. For the CXO spectrum 
(Epoch 1, Fig.[l]3l), which has fewer counts, we derived 
our model fit with two components as above, but without 
a Fel line (as none could be seen in the spectrum), and 
we fixed the iron abundance to the XMM-model value. 

For Cl Hyi and EY Hya, we fit the EPIC spectra with a 
1-T APEC model (Fig.[]Ji2,a3), allowing [Fe/H] to vary. 
A gaussian Fe I line component was also included for fit¬ 
ting EYHya, in the same manner as for the modeling 
of the Y Gem XMM spectrum. The derived [Fe/H] val¬ 
ues are not significantly different from the solar value. 
For the lower S/N ACIS spectrum of EYHya (Fig.[T|32), 
we fitted 1-T models fixing the value of [Fe/H] to be 
the same as in the EPIC fit, and derived the required 
Tx and NH. For the even lower S/N ACIS spectrum of 
CIHyi (Fig.[l]33), we fixed both Tx and [Fe/H] to be the 
same as in the EPIC fit, and derived the required NH. 
Our modeling shows significant variations in NH for each 
source. 

These data, together with previous studies (Ramstedt 
et al. 2012 [RMKV12], Kastner & Soker 2004 [KS04]), 
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show that X-ray emission is found only in fuvAGB stars - 
these include 3 objects from our study (CIHyi, EYHya, 
Y Gem) and 2 from archival data, R Uma and T Dra 
(RMKV12). We exclude the fuvAGB star SSVir from 
this list, since its tentative detection, 4.8 db 1.3 counts 
ks -1 , may be an artifact of optical loading (RMKV12). 
For the set of AGB stars not detected in X-ray emission 
but with adequate exposure times, >0.5 ks (thus allowing 
meaningful upper limits) and that were observed with 
GALEX in the FUV, we have 2 objects from archival 
data, RTEri and TXCam (RMKV12, KS04) in addition 
to the 3 non-detected objects from our study (NUPav, 
delOl Aps, and V Hya). However, as TX Cam and RT Eri 
are not fuvAGB stars (we find 3cr upper limits of 8.5 and 
12.4 fi Jy, respectively, for their FUV fluxes), these are 
not considered further. 

Amongst the well-known symbiotic systems with AGB 
primaries, Mira, RAqr and CHCyg, only Mira has been 
observed in the FUV with GALEX, is both an X-ray 
emitter (Karovska et al. 2005) and a fuvAGB star (Se- 
talll), and believed to have a WD companion (Sokoloski 
& Bildsten, 2010). 

In Fig. EH, we plot the unabsorbed, average, X-ray flux 
versus the observed average (GALEX) FUV flux for fu¬ 
vAGB stars (and Mira) - it appears that the FUV flux 
is not a good indicator of the X-ray emission. For the 
FUV data, the error bars show the flux ranges when 
data are available for more than one epoch. For the 
X-ray data, the error bars are as follows : (i) for the 
3 detected sources in this paper, the errors bars show 
the range covered by the X-ray fluxes over two epochs, 
(ii) for the 2 detected sources from RMKV12, the error 
bars show the range in the unabsorbed flux assuming 
log(NH) ~ 21.5 =b 0.4 (as NH is not well constrained 
in their models, we assume a typical NH value and a 
conservatively large uncertainty), (iii) for each of the 3 
non-detected sources in this paper, we use the Lx ver¬ 
sus NH relationship in Fig. 4 of RMKV12 for T Dra, and 
log(NH) rsj 21.5 zb 0.4, to convert the upper limit on the 
X-ray count rate to an Lx upper limit. The FUV fluxes 
have not been corrected for dust extinction (intrinsic or 
extrinsic to the source) as we do not have a reliable way 
of estimating these - however, it is unlikely that account¬ 
ing for these would affect our conclusion above. Since our 
co-eval observations show at least a modest correlation 
between the variations of the X-ray and UV fluxes, im¬ 
plying a close relationship between the X-ray and FUV 
emission, the observed absence of a correlation in non¬ 
coeval observations is likely due to the fact that both are 
variable. 

In Fig.Eb, we plot the average X-ray (i.e., 0.3-10 kev) 
luminosity versus the GALEX FUV/NUV flux ratio, 
Rfuv/nuv For stars with multi-epoch UV observations, 
the error bars cover the full range between minimum and 
maximum values of Rfuv/nuv We find Rf uv /nuv >0.17 
in the 5 fuvAGB stars with X-ray emission (and Mira.) 
Amongst the 3 fuvAGB stars without X-ray emission, 
Rfuv/nuv < 0-12 for NUPav and delOlAps; for VHya, 
Rfuv/nuv = 1.1. Thus Rf uv /nuv appears to be a better 
indicator of X-ray emission than the FUV flux - but we 
must caution that this is a tentative result since the total 
number of fuvAGB sources in this plot is small. 

3. DISCUSSION 


We now discuss several models for the X-ray emission 
that we have observed from fuvAGB stars - all but one 
of these require binarity (note that at present there is no 
other evidence of binarity in our targets). 

3.1. X-ray and UV Variability in fuvAGB Stars 

The X-ray and UV variability provide useful con¬ 
straints on the nature and diagnostics of the emitting 
region. The observed variability is perhaps the strongest 
indication that these emissions are related to the pres¬ 
ence of a binary companion. 

Although the short-term variations in EYHya and 
CIHyi don’t show obvious periodicity, the variability 
time-scales are similar to that in Y Gem, hence it is plau¬ 
sible that the variability mechanism is the same in all 
three objects. 

It is unlikely that the X-ray variability is due to flare 
activity since the fuvAGB stars’ X-ray light curves do 
not show the temporal behaviour seen in YSO flares 
that have rapid rise times followed by slower exponential 
decay. Secondly, in flares the hardness ratio increases 
during the flare peaks, but for Y Gem, where the XMM 
count rate is high enough to determine the hardness- 
ratio variations, we find an anti-correlation between these 
two quantities (Fig. [5]). We defined the hardness ratio as 
(.H — S)/(H + S'), where H and S are the fluxes in the 
energy bands 4-10 keV and 2-4 keV. We excluded the 
flux at energies shorter than 2keV in defining this ra¬ 
tio because the former has much larger errors due to the 
smaller number of counts. 

The hardness ratio variations in Y Gem are consistent 
with the X-ray changes occuring predominantly as a re¬ 
sult of changes in NH (for EY Hya and Cl Hya, the count 
rates are too low for such hardness-ratio light-curve anal¬ 
ysis.) We note that our modeling of the X-ray spectra 
shows that NH varies significantly on long time scales in 
all 3 X-ray fuvAGB stars. It is thus likely that NH vari¬ 
ations are the dominant cause of short and long-term 
X-ray variations in fuvAGB stars, which may be pro¬ 
duced by changes in viewing geometry (e.g., due to a 
warped rotating accretion disk), the presence of variable 
accretion streams, or both. 

3.2. Models for X-ray Emission from fuvAGB Stars 

3.2.1. Model 1. Accretion Disk around a Companion 

If we assume that the periodicity in Y Gem is asso¬ 
ciated with the orbital period of a central binary, then 
using Kepler’s law, we find that the semi-major-axis is, 
a<0.0036 — 0.0046 AU or (5.4 — 6.8) x 10 10 cm, assum¬ 
ing primary and companion masses of M p = 1 — 3 M 0 , 
M c < 1 Mq. Since such values for a would place the com¬ 
panion deep inside the AGB star’s photosphere, the X- 
ray variation is unlikely to be associated with the orbital 
period of the binary. 

The ^1.3 hr time-scale is similar to the period of ma¬ 
terial orbiting close to the inner radius of an accre¬ 
tion disk around a sub-solar mass companion, i.e., with 
M c < 0.3 5 Mq (implying a< 3 x 10 10 cm); for larger masses, 
the orbit radius becomes smaller than the stellar radius. 
Y Gem’s X-ray and UV variability therefore suggests that 
these emissions arise at or near the magnetospheric ra¬ 
dius in a truncated disk, or the boundary layer between 
the disk and star. In either case, the relatively large emis- 
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sion measures derived from our APEC modelling imply 
large densities, n, in the emitting region. Assuming that 
the emission comes from a torus of inner radius r\ n and 
circular cross-section with fractional radius Sr, we find 
that 

n = 2 x 10 11 cm- 3 (EM/10 54 cm" 3 ) 1 / 2 

x (3 x 10 10 cm/r in ) 3/2 (0.2/S r) 1/2 . 1 

The magnetic field required to confine such a dense, hot 
plasma is correspondingly large since the magnetic en¬ 
ergy density must be greater than or equal to the ther¬ 
mal energy density. We find that the magnetic field at 
the inner radius, is, 

B in > 260 G (n/2 x 10 11 cm- 3 )(Tx/10 8 K). (2) 

For a truncated disk, the magnetospheric radius, r ma g , 
depends on the magnetic field of the companion, B c , and 
the accretion rate ( M acc ) (see Eqn. 8 of Stute & Ca- 
menzind 2005). Assuming the same geometry for the 
X-ray emitting region as above, and the derived value of 
Bi n (260 G), extrapolated to the stellar surface assum¬ 
ing a dipole field, we find that, r ma ^ = 2.9 x 10 10 cm, if 
M acc ~ 1O- 9 M 0 yr -1 and M c = 0.25M 0 . The mag¬ 
netic field at the stellar surface is, £>* = 560 G. The 
derived value of r map is relatively insensitive to the un¬ 
known accretion rate, since r ma ^ oc M a cP • 

Since the values of r mag and r^ n are so similar, then 
given the uncertainties in the estimates of both, due to 
our assumptions on the geometry of the emitting region, 
the magnetic field geometry, and the accretion rate, it is 
possible that the disk is not truncated, but is connected 
to the star via a boundary layer. 

3.2.2. Model 2. Accretion Onto a Companion 

The generally accepted model for accretion from the 
disk into YSOs (the magnetospheric model) envisions 
plasma being channeled into magnetic flux tubes and 
ramming onto the star at essentially free-fall speed (e.g., 
Favata 2005). The X-ray temperature of the shocked 
gas resulting from such a process is, T s = 3.44 x 
10 6 K(M C /R c ) (from Eqn. 9 of Calvet & Gullbring 
1998), where the companion mass and radius, M c and R c , 
are in solar units. Thus, for a Sun-like companion, T acc is 
significantly lower than observed; for more massive dwarf 
companions, even upto B5, M c /R c < 1.5. Consider¬ 
ing substellar companions, only the most massive brown- 
dwarfs (BDs) provide a somewhat larger value of M c /R c , 
e.g., for a BD with mass 0.1 M 0 , M c /R c = 1.5. Only for 
white-dwarf (WD) companions, with M c /R c > 10 — 45, 
the observed values of Tx (~ 35 — 160 K) can be easily 
produced as a result of accretion. 

However, none of our X-ray emitting AGB stars is 
known to be a symbiotic star with a WD companion. For 
EY Hya, an optical spectrum taken with the Palomar 5m 
telescope shows no emission lines in the 3889-5436 A re¬ 
gion (K. Findeisen, private communication). But we 
cannot rule out the possibility that fuvAGB stars have 
WD companions that are not hot enough to ionize a de¬ 
tectable amount of gas and produce emission lines (e.g., 
the Balmer lines or forbidden lines such as [OI] A 6300 or 
[Nil] AA 6549, 6583). If so, these WD companions must 
have cooled to effective temperatures T e ff < 20, 000 K. 


However, in the case of Y Gem, Setalll have shown that 
its intense UV emission requires a much larger surface 
area than that of a WD, hence its UV emission cannot 
be due to accretion onto a WD. 

However, if an adequate mass of accreting gas can be 
pre-accelerated to high speeds (v s ) by magnetic reconnec¬ 
tion (e.g., de Gouveia Dal Pino et al. 2014, Hamaguchi 
et al. 2012), then the shock temperature is, using Eqn. 9 
of Calvet & Gullbring (1998), T s = (3/16 ks) limner]. ~ 
70 MK (f s /1500 km s -1 ) 2 , assuming a mean molecular 
weight, /i = 1.33, and is not constrained by the mass-to- 
radius ratio of the companion. 

Hamaguchi et al. (2012) propose that such “accretion- 
induced magnetic reconnection” generates the X-ray 
emission in the YSO V1647 Ori; the differential rotation 
between star and its accretion disk shears the star’s mag¬ 
netic field, causing the field lines to twist and continu¬ 
ously reconnect. 

The quasi-periodic variability in the X-ray light curve 
of V1647 Ori, rather similar to what we find for the fu¬ 
vAGB stars, is a by-product of this mechanism due to 
the production of hot-spots on the rotating central star, 
where, matter that is accelerated to high speeds by the 
magnetic reconnection collides with the star to produces 
X-ray emission with high Tx. It is possible that such 
a model applies to Y Gem (and possibly EY Hya and 
CIHyi) as well. 

3.2.3. Model 3: Stellar Coronal Emission 

X-ray emission from stellar coronae has been exten¬ 
sively studied (e.g., Giidel 2004). The relatively high 
values of Tx that we have found in our fuvAGB sources 
argues against the X-ray emission coming from such coro¬ 
nae, which typically show values of Tx in the range 
~ (2 — 10) MK and rarely as high as the lowest Tx value 
in our sample, ~ 40MK (e.g., Schmitt et al. 1990). We 
also note that stellar coronae are usually accompanied 
by stellar chromospheres, and the latter produce strong 
Ca IIH & K emission lines - however the optical spectrum 
of EYHya (K. Findeisen, private communication) shows 
no such emission. We discuss coronal models in more 
detail below, and present additional arguments against 
them as being the likely source of X-ray emission in fu¬ 
vAGB stars. 

(a) From a main-sequence companion — “Saturated” coro¬ 
nae around low-mass (< 1M 0 ) main sequence stars 
reach a maximum of Lx/Lc^ 10 -3 (Pizzolato et al. 
2003). More massive stars (1.1-1.29M 0 ) follow a lower 
ratio, Lx/Lc^ 10 -3,9 . Our observed values of Lx he 
in the range 0.002 — O.2L 0 . Then, if this X-ray lu¬ 
minosity comes from coronal emission from a main se¬ 
quence companion, the implied companion luminosity is 
Lc> (2 — 50) L 0 . But since Lc> 2 L 0 , the companions 
must have masses >1.1 M 0 . Thus the implied compan¬ 
ion luminosities are even higher, Lc> 15-500 L 0 , implying 
masses Mc>12M 0 , which is not allowed, since the com¬ 
panion mass must be less than that of the primary AGB 
star (1-8 M 0 ). 

(b) From the AGB star — Since plasma at temperatures 
> 0.5 keV cannot be gravitationally confined on the sur¬ 
face of giants and supergiants (Rosner, Golub & Vaiana 
1985), the very high values of Tx that we have found 
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imply that the emitting plasma in fuvAGB stars, must 
be magnetically confined (and heated). A universal rela¬ 
tionship between magnetic flux and the power dissipated 
through coronal heating is suggested by the linear de¬ 
pendence between the total unsigned magnetic flux, 
and Lx for the Sun and active dwarf stars (Pevtsov et 
al. 2003) over 12 orders of magnitude. Applying this 
relationship to our fuvAGB stars, with Lx = (0.92 — 
42) x 10 31 ergs -1 , we get <pB = (0.92 — 42) x 10 27 G cm 2 , 
implying, for a AGB star with radius 1 AU, an average 
magnetic field B av = (0.3 - 15) G / f ag b, where f ag b is 
the fraction of the AGB’s star’s surface covered by the 
B field (fagb ^ 1 for a large-scale field and f ag b « 1 
for localized fields). Since we also need to produce sub¬ 
stantial amounts of plasma implied by our derived val¬ 
ues of the X-ray emission measure, around the fuvAGB 
stars, and since these objects are cool, late-type stars, 
with three of them known to drive cool, dusty molecular 
winds (EYHya, TDra, RUma), the surface filling factor 
of the coronal gas (and thus f ag b) is likely to be relatively 
small. 

The linear relationship between Lx-0# appears to satu¬ 
rate at the largest values of Lx (i.e., at Lx 10 30 ergs 1 , 
in a region populated by T Tauri stars) with 5/6 objects 
showing Lx values a factor of 10 or more below the lin¬ 
ear fit. Hence the value of B av implied by Lx in fuvAGB 
stars may be higher by a factor 10 or more. Recent stud¬ 
ies show that such values of B av are plausible: e.g., mag¬ 
netic fields strengths of 0.3 — 6.9 (15.8 — 1945) G have 
been inferred for the stellar surfaces in 3 AGB stars by 
extrapolating the field-strength from observations of po¬ 
larization in the 22 GHz H 2 O maser line, assuming an r -1 
(r -3 ) radial variation, appropriate for a toroidal (dipole) 
field (Leal-Ferreira et al. 2013). 

However, a sensitive XMM search for X-ray emission in 
two AGB stars with known or suspected strong B-fields 
yielded null detections (KS04), making coronal emission 
from the AGB star a less likely candidate for producing 
the observed X-ray emission. 


4. CONCLUDING REMARKS 

In a survey for X-ray emission from 6 AGB stars with 
far-ultraviolet excesses, we detected 3/6 objects with rel¬ 
atively high X-ray luminosities and plasma temperatures; 
furthermore the X-ray flux in each of these is variable 
on hour-long time-scales. In contrast, a similar sensitive 
search for X-rays in 2 AGB stars with evidence for mag¬ 
netic fields, was unsuccessful, and these stars show no 
FUV emission. We believe therefore that the most likely 
scenario for the origin of the X-ray (and FUV) emission 
involves accretion activity around a companion star, with 
confinement by strong magnetic fields associated with 
the companion and/or an accretion disk around it. How¬ 
ever, an extended survey of AGB stars with and without 
UV excesses (including those with evidence of magnetic 
fields) is needed to help us robustly test this hypothesis. 
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TABLE 1 
Observations log 


XMM-Newton observations 


Target 

Camera 

Filter 

Eff. Exp 

V 

Date 

X-Ray 

Spect. 




Time (sec) 

(mag) 

dd/mm/yyyy 

Emiss.? 

Type 

Cl Hyi 

pn 

Medium 

7897 

9.3 

03/10/2013 

Yes 

M6 


MOS1 

Medium 

13843 






MOS2 

Medium 

14006 






OM 

UVM2 

4000 






OM 

UVW2 

11700 





Y Gem 

pn 

Medium 

6238 

9.09 

31/03/2014 

Yes 

M8 


MOS1 

Medium 

8250 






MOS2 

Medium 

8510 






OM 

UVM2 

4000 






OM 

UVW2 

4000 





EY Hya 

pn 

Medium 

16839 

9.37 

06/11/2013 

Yes 

M7 


MOS1 

Medium 

17284 






MOS2 

Medium 

17968 






OM 

UVM2 

4400 






OM 

UVW2 

12800 





delOl Aps 

pn 

Thick 

11539 

4.76 

18-19/08/2013 

No 

M5 


MOS1 

Thick 

13166 






MOS2 

Thick 

13172 





NU Pav 

pn 

Thick 

10039 

5.06 

26-27/09/2013 

No 

M6 


MOS1 

Thick 

11701 






MOS2 

Thick 

11672 





Chandra observations 

Target 

Camera 


Eff. Exp 

V 

Date 

X-Ray 

Spect. 




Time (sec) 

(mag) 

dd/mm/yyyy 

Emiss.? 

Type 

Cl Hyi 

ACIS-S 


9838 

9.3 

24/07/2014 

Yes 

M6 

Y Gem 

ACIS-S 


10335 

9.09 

15/12/2013 

Yes 

M8 

EY Hya 

ACIS-S 


9839 

9.37 

25/12/2013 

Yes 

M7 

V Hya 

ACIS-S 


9660 

9.70 

18/12/2013 

No 

C 







TABLE 2 

X-Ray Emission Properties of fuvAGB Stars 


Detected Sources: Models 


Target 

NH 

log(Tx) 

log(EM) 

Fx 1 

Lx 2 

D 3 

flux[FeI(6.4)] 4 

[Fe/H] 5 

x 26 


10 22 cm -2 

K 

cm -3 

erg cm -2 s -1 

io- 3 l @ 

kpc 

phot s -1 cm -2 

(dex) 


CI Hyi/XMM 

0.33 

7.87 

54.2 

7.47 X 10" 13 

7.8 

0.58 


0.12 

1.2 

range 7 

0.05, -0.04 

0.04, -0.07 

0.03, -0.03 

±8.3 x 10- 16 




0.13,-0.17 


CIHyi/CXO 

2.8 

7 . 87 s 

54.2 

4.06 x 10“ 13 

7.9 



0.12 s 

1.2 

range 7 

0.76, -0.59 


0.08, -0.09 

±3.5 x 10“ 15 






EY Hya/XMM 

0.095 

7.57 

53.7 

5.56 x 10" 13 

2.0 

0.35 

6.7 x 10“ 7 

-0.16 

1.35 

range 7 

0.009, -0.009 

0.05, -0.04 

0.02, -0.02 

±3.2 x 10“ 16 



±6 x 10" 7 

0.11,-0.13 


EY Hya/CXO 

0.05 

7.74 

53.7 

5.85 x 10- 13 

2.4 



- 0 . 16 s 

0.70 

range 7 

0.04, -0.04 

0.13,-0.11 

0.04, -0.04 

±5.8 x 10“ 15 






Y Gem/XMM 

8.13 9 

8.1 

55.6 9 

1.05 x lO -11 

226 

0.58 


-0.10 

1.1 

range 7 

0.40, -0.38 

0.03, -0.04 

0.02, -0.02 

±1.6 x 10“ 15 




0.06, -0.07 



0.043 10 

8 . 1 s 

53.7 10 




2.4 x 10“ 5 



range 7 

0.02, -0.015 


0.05, -0.05 




±4.8 x 10“ 6 



Y Gem/CXO 

15.6 9 

8.2 

55.2 9 

3.1 x 10- 12 

115 



- 0 . 10 s 

1.2 

range 7 

16.6, -5.8 

0.52,-0.13 

0.14,-0.21 

±7.5 x IQ" 14 







3.04 10 

8 ,.£ 8 

54.6 10 







range 7 

1.8, -1.5 


0.35, -0.49 







Non-Detected Sources 

delOl Aps/XMM 





< 0.039 11 

0.23 




NU Pav/XMM 





< 0.030 11 

0.16 




V Hya/CXO 





< 0.090 11 

0.40 





a The observed X-ray flux in the 0.03-10 kev range 
b The intrinsic X-ray luminosity in the 0.03-10 kev range 

c Distanees: from Hipparcos parallax; if unavailable, using MK=-7.6 for late-M semi-regular stars, as in Kahane & Jura (1994) 

d FeI 6.4 kev line flux (a Gaussian line-shape was assumed and the line energy and width (FWHM) were fixed at 6.4 keV and 0.002 keV during fitting) 
logarithm of the ratio of star’s Fe abundance to that of the Sun 
f reduced y 2 value 

g ±l (7,-1 cr values for NH, log(Tx), log(EM), Fx, and the Fel(6.4) line flux 
^Parameter value (in italics) was fixed, hence no uncertainties are provided 
i NH(l), EM(2) 

JNH(2), EM(2) 
k 3cr upper limit 


-A 



00 



Energy [keV] 


Energy [keV] 


Energy [keV] 


Fig. 1 . — X-ray spectra (colored curves) and model fits (black curves) of the fuvAGB stars Y Gem, EYHya and CIHyi. Panels rows show XMM/EPIC (pn: red, 
MOS1: green, MOS2:blue), and CXO (ACIS-S). Spectra shown include the background; the ISIS model fitting procedure takes the latter into account (see text for details). 
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EYHya 


(al) (a2) 



CIHyi 


(a3) 



17 17.5 18 18.5 19 19.5 


Time (hrs since JD 2456862) 


Fig. 2.— X-ray and UV light-curves of the fuvAGB stars YGem, EYHya and CIHyi. The energy range used to extract the X-ray curves was 0.3-10keV for XMM, 
and 0.3-7keV for CXO. Panels show data from XMM (EPIC=pn+MOSl+MOS2: red, MOS=MOSl+MOS2: green, UVM2: black squares, UVW2: black circles), and CXO 
(ACIS-S). The EPIC, MOS and UVW2 data have been respectively re-scaled as follows: 80,130, 2.5 (YGem), 2, 5, 0.8 (EYHya), 20, 30, 2.5 (CIHyi). A sinusoidal fit (by eye) 
with period, P=1.35hr is shown for Y Gem. All data are background-subtracted; error bars are Tier. In order to facilitate the comparison between the XMM and CXO 
light curves, the widths of the two panels for each source have been adjusted so that similar time-intervals have similar lengths (on the horzontal axes). 


co 











Counts/s/keV 


10 
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Energy [keV] 

Fig. 3.— The Fe line region in the Y Gem EPIC/pn spectrum. The expected locations of the Fel line and the line complexes of the 
coronal FeXXV and FeXXVI lines, are marked. 
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(b) 

100 



CO 

I 

o 


co 

O 

c 


E 




cc 

I 

X 


10 


i 


o.i 


o.oi 


Y Gem 


o 

T Dra 


R Uma 

G-© 

Cl Hyi 


G-G 

EY Hya 

Mira-flare 

o 

o- ©Mira 


delOl Aps 
NU Pav 


G-G 

V Hya 


0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

FUV Flux (gjy) / NUV flux (gjy) 


Fig. 4. — (a) The average (unabsorbed) X-ray flux versus the average (absorbed) FUV flux for fuvAGB stars (green symbols) and Mira, 
a symbiotic star with an AGB primary (blue symbol). Upper limits (3cr) are shown for fuvagb stars not detected in X-rays (red open 
symbols), (b) Average (unabsorbed) X-ray luminosity versus the FUV-to-NUV flux ratio, Rf uv / n uv> f° r fuvAGB stars, with green & blue 
symbols representing X-ray detections (red symbols show 3cr upper limits). In the case of the symbiotic star, Mira, the plot shows its 
quiescent X-ray luminosity versus the Rf uv /nuv range, as well as the X-ray luminosity during a flare on 2003, Dec 6 versus the Rf U v/nuv 
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Fig. 5. — Light curves of the X-ray EPIC/MOS count-rate (green) and the hardness ratio (black), defined as (H — S)/(H + S), where H 
and S are the fluxes in the energy bands 4-10 keV and 2-4 keV, for Y Gem (error bars are d=l<r.) 


























